The Dark Energy Cosmic Clock: A New Way to Parametrise the Equation of State 
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We propose a completely new parametrisation of the dark energy equation of state, which uses the 
dark energy density, Q,. as a cosmic clock. We expand the equation of state in a series of orthogonal 
polynomials, with fie as the expansion parameter and determine the expansion coefficients by fitting 
to SNIa and H{z) data. Assuming that f2c is a monotonic function of time, we show that our 
parametrisation performs better than the popular Chevallier-Polarski-Linder (CPL) and Gerke and 
Efstathiou (GE) parametrisations, and we demonstrate that it is robust to the choice of prior. 
Expanding in orthogonal polynomials allows us to relate models of dark energy directly to our 
parametrisation, which we illustrate by placing constraints on the expansion coefficients extracted 
from two popular quintessence models. Finally, we comment on how this parametrisation could be 
modified to accommodate high redshift data, where any non-monotonicity of flc would need to be 
accounted for. 



I. INTRODUCTION 

During the last decade a vast amount of cosniological 
data has been collected, which indicate that a mysterious 
form of dark energy is driving an accelerated expansion 
of the universe [11-0] . The simplest explanation for dark 
energy is Einstein's cosmological constant. A, which has 
a constant equation of state We = Pel Pe — —1- How- 
ever, A suffers from several short-comings, including the 
fine-tuning and coincidence problems (see e.g. [5] for a 
review) . 

These thorny issues surrounding the cosmological con- 
stant have prompted investigation into alternative mod- 
els such as quintessence [SHE], Chaplygin gas |9j, modi- 
fied gravity [TU] , holographic dark energy [IT] , and many 
others, all which promote A to a dynamical degree of 
freedom with a time-varying effective equation of state. 
For a review of models of dynamical dark energy see [H] . 
Instead of appealing to a fundamental theory to describe 
We{t), we can attempt to reconstruct its properties in a 
model independent way by proposing a functional form 
for Wc{t) and fit this directly to observation. Constraints 
may be placed on the dark energy equation of state once 
a parametrisation has been adopted. For example, by as- 
suming a constant equation of state, We = w, the authors 
of [m report w = -1.008 ± 0.085 at 68% CL, consistent 
with the ACDM concordance cosmology. Very recently, 
via a principle components approach, the authors of |14| 
find that Baryon Acoustic Oscillation (BAO) data sug- 
gest deviations w{z) < — 1 above one standard deviation 
at redshifts z ~ 0.25. 

To the best of our knowledge, all of the Wc{t) 
parametrisations that have been proposed in the liter- 
ature to date, use redshift z (or equivalently scale factor 
a) as the 'time' variable, i.e., wdt) — f{z). For example. 
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the Chevallier-Polarski-Linder (CPL) parametrization, 
first discussed in [T5] and reintroduced in [IB] uses a poly- 
nomial fitting function in redshift space, whilst in |17j a 
logarithmic expansion in z was proposed. Motived by the 
dynamics of quintessence models, Ref. [TBj introduced a 
parametrisation dependent on five parameters which is 
able to reproduce the time evolution of Wc across a wide 
range of redshifts for a variety of different quintessence 
models. Other studies have analysed our ability to re- 
construct the quintessence potential V{(j)) [HHUj- 

In this paper, we introduce a completely new 
parametrisation of We{t), which uses the dimensionless 
dark energy density fraction ric(0 — Pc{t)/'iH{t)^ as a 
cosmic clock. The idea is to expand wd^c) in orthogonal 
polynomials, with flc as the expansion parameter: 
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This is in similar spirit to [55] where the authors were 
interested in parametrising the evolution of small scale 
density perturbations. Such an expansion has several 
advantages. Perhaps most importantly, Q^ is a physical 
quantity, directly related to the properties of dark en- 
ergy. Furthermore, assuming < f7o < 1, it makes for an 
ideal expansion parameter since ilo is a naturally small 
number. Expanding Wc in terms of orthogonal polyno- 
mials has been carried out before, see for example |23) . 
but with redshift as the expansion parameter. 

So long as ilo(i) remains a monotonic function over 
the epoch of interest it may be used as a perfectly good 
cosmic clock. As noted in [24] it is natural to consider 
that flo increases monotonically through most of cosmic 
history, an assumption that is well motivated by various 
astrophysical constraints: 

• Big Bang Nucleosynthesis (BBN): Vtc < 0.045 at 
z - lO^o [25^ , 

• Galaxy formation epoch: J7o < 0.5 at z ^ 2 — 4 , 

• Present day: ^c ^ 0.68 with w^ ^ —1 at z = '26j . 



This monotonicity was recognised in [57] and [55], 
where the authors proposed a dark energy parametrisa- 
tion in terms of three parameters: the dark energy equa- 
tion of state today wq, the amount of dark energy today 
fi^j Eind the amount of dark energy at early times il^ to 
which it asymptotes at high redshift. Our parametrisa- 
tion is similar to this in the sense that it directly relates 
Wc{t) to ilc(Oi however our parametrisation does not rely 
on scale factor/redshift time, unlike Refs. [571155] . 

We also know that the transition from matter to dark 
energy domination (and hence the start of cosmic acceler- 
ation), occurred at low redshift z sa 1 — 2. This is a region 
that surveys such as the Dark Energy Survey (DES) will 
be sensitive to, and will depend crucially on fie where the 
transition occurs. Our parametrisation is ideally suited 
to such a case since we are expanding in a small param- 
eter around the time of transition. We can in principle 
use our parametrisation to reconstruct the equation of 
state around that region and in doing so provide another 
handle on understanding the nature of dynamical dark 
energy. 




II. THE DARK ENERGY COSMIC CLOCK 

Assuming a spatially flat Friedman-Robertson- Walker 
(FRW) metric, the background equations for any theory 
of gravity can be recast in the usual form as used in GR. 
The Friedman equation reads 



H^t) 




(2) 



where a and H are the scale factor and Hubble func- 
tion respectively, pc is the dark energy density (which 
may in general be a function of additional degrees of 
freedom), and pi are the energy densities of the other 
possible components, including matter pm and radiation 
Pj. An overdot denotes differentiation with respect to 
cosmic time t, and we use natural units throughout, with 
SttG = TOpj^ — 1. Regardless of the theory of gravity, one 
may always treat dark energy as a standard fluid with a 
time-varying equation of state Wc{t) — Pe(t) / pe(t) , sub- 
ject to energy conservation p^ + 3iJ(l + 'Wc)pc = 0. Once 
the theory of gravity and the other components pi are 
specified, Wc may be subject to additional equations de- 
scribing its time evolution. 

The energy conservation equation for pc can be easily 
rewritten as 



n. 
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(3) 



where the sum is taken over all other cosmic compo- 
nents i. In a universe containing dark energy, pressure- 
less matter {wn 

whenever Wc(l 



= 0), and radiation {w^ — 3); ^^ have 
1, and there are turning points in fl^ 
flc) = 5^7 (assuming that _ff > at 



FIG. 1. The exact equation of state for four quintessence mod- 
els, li)thoory(fio) = Pc/ pc, (solid lines) and the corresponding 
reconstructed equation of state, (dashed lines) obtained by 
building an interpolating polynomial Wo(^g) = X^„ WnUn{^e) 
over the interval fJe = ^c{z = 1.75), ilo], keeping the first 
three terms in the expansion. Time advances from left to 
right, and the vertical dashed line denotes f2o, the value of ilo 
at 2 = 0. 



all times). Hence, monotonicity of ilc is broken during 
the radiation dominated era whenever Wc — ^ and in 
the matter era whenever w^ = 0. At these points, the 
dark energy clock would 'stop ticking', and our expan- 
sion Eq. (fTl) would break down. 

Many scalar field dark energy (quintessence) models 
possess scaling solutions on which the scalar field energy 
density tracks that of the dominant background fluid, Qi. 
For example, in the case of a single exponential scalar 
field potential, V{(l>) = Vbe"^"^, scaling solutions exist 
whenever A^ > 3{wi + 1). On these solutions the scalar 
field dark energy equation of state mimics the evolution 
of the dominant background fiuid, Wc = Wi [291 130]. This 
is an explicit example of where we can expect the dark 
energy clock to break down. 

In this paper, we use our dark energy clock parametri- 
sation to fit to low redshift background expansion data 
only, and so we neglect radiation. The only fixed point of 
Eq. ([3]) that concerns us then is the matter era Wc = so- 
lution. Hence, for a given set of expansion coefficients w„, 
ric will possess turning points whenever '^^ i(;„P„(i7o) = 
0. We will return to this important point in the next sec- 
tion. Our choice of P„ is the Chebyshev polynomials of 
the second kind. By defining a suitable inner product we 
may choose the interval, fi^ = [fi™'", 17™'''^], over which 
the polynomials are orthogonal. We denote these shifted 



Chebyshev polynomials of the second kind by C/„(17c), 
and write 



,(rio) = y^Wn^«(^c 



(4) 



The n zeros (nodes) of the C/„ are useful for interpolation 
because the resulting interpolation polynomial minimizes 



Runge's phenomenon (the problem of oscillation of the 
interpolating polynomial near to the edges of the inter- 
val). The properties of the C/„ and how they are related 
to the standard Chebyshev polynomials of the second 
kind [/„, are given in Appendix \K\ 



From the C/„ orthogonality condition, Eq. ( A4), we can 



extract the expansion coefficients Wn , given any smoothly 
varying monotonic w{ne): 



w„ 
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Hence, any dark energy model which predicts a mono- 
tonic equation of state i(;thoory(^c)j can be directly re- 
lated to our parametrisation. In Fig. IT] we show the equa- 
tion of state, u;thoory(^c), for four different quintessence 
potentials V{4)): 



lEXP 


Foe^* [291130], 


2EXP 


yo[e^"^ + e^^'^] [SH 
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Vo [Ai + (A2 - (bf 
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Vo^+V^^e^'/^ [31] . 
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In the same figure, we also plot w^^^c), the 'recon- 
structed' equation of state built from shifted Chebyshev 
polynomials orthogonal on the interval Qo = [^c{z = 
1.75), ilo], where we keep the first three terms in the 
expansion, and Jig = ^o{z = 0). The height of the po- 
tential (set by Vo) is adjusted to give ilc — 0.72 for each 
model. Since the dynamics of the scalar field is different 
for each model, the evolution of the dark energy density 
is also different, and so the lower limit ^e{z = 1.75), will 
in general be different from model to modeljj 

As can be seen from the figure, our expansion does 
well in capturing the evolution of w^. The rapid 
oscillations seen for the AS model are induced by the 
late time oscillations of the field about its minimum. 
Our parametrisation is unable to resolve the individual 
oscillations (to do so would require retaining a large 
number of terms in the expansion), but it does capture 
the average behaviour rather well. Furthermore, for 
all four models, the series rapidly converges. For 
example for the lEXP model we find: Wq — —0.876, 
wi — 0.050, W2 — 0.005. In general of course, the rate 
of convergence will depend upon the behaviour of the 
function We(i^c) that we are trying to reconstruct. If, 
as various astrophysical constraints suggest, Qc becomes 
less important at high redshift, our expansion should in 



^ A cosmological constant (wc = —1), maximally decreases 
Cc(z) with increasing redshift relative to the non— phantom 
quintessence models considered here. For ACDM, with CaC^ = 
0) = 0.7, we have Qa=0.1 at 2 ~ 1.75. 



principle always converge, as higher order terms f^JJ 
large n will become negligible. 



for 



Having motived our Wc parametrisation, we now turn 
our attention to assessing its performance by fitting to 
background expansion data. To aid this analysis we first 
present analytic solutions to the background equations 
of motion for dark energy which are valid up to second 
order in the expansion of We(f^e)- 

Scale factor, a(flc)'. We begin with our equation of 
state parametrisation in terms of Chebyshev polynomi- 
als, Eq. Q, which may be rewritten at second order as: 



Wei^e) = Wo + Wifle + W2^l 



(6) 



The coefficients, wq, wi and 'W2 are combinations of the 
Wn and are given in Appendix |A] Neglecting radiation, 
Eq. PI may be rearranged to give: 



"da 



dn. 



r2o(l — ^c)(wo + WiQc + 'W2^c) 



(7) 

where the lower limits of the integrals (subscript 0) de- 
note the value today, and we have written fleiz = 0) = 
flo . We remind the reader that wq is not the value of w^ 
today, but is the zeroth order expansion parameter. The 
LHS of Eq. (ItI) is trivial, and the RHS may be expanded 
using partial fractions and the resulting terms integrated 
separately. We find: 
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,F(n„ 



(8) 

where we have set ag — 1. The powers a, /3 and 7 are 
combinations of the w„, and correspond to coefficients of 
the partial fraction expansion of Eq. (ItI). As a result, if 
any one or more of the Wc{^c) expansion coefficients wq, 
Wi or W2 in Eq. nh is exactly zero, then the powers a, 
/? and 7 and the function F{nc) will change. In the case 



where none of the w„ are zero: 

1 + Wq ^ 1 + Wt 

wt 

, Wt — WQ + W1+W2, (9) 
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, /3 

1 Wi + W2 

2u'o Wt 



and 



F(r2o) = [wi{wi + W2) - 2W0W2] 



X arctan 



2w2q{^o - ^c) 



(2w2f^o + 'Wi){2w2^c + wi) + g2 



(10) 



where q = •\/4u;oW2 — w\. The solutions corresponding 
to the six different cases of zero Wq, Wi or W2 are 
hstcd in Table |IV] of Appendix [B] Notice that the 
solution (Is]) breaks down at flc — 0, ^c — ^ and 
Wo(^c) = Wo + "u^i^c + ''i'2^c =0- These are fixed points 
of Eq. (3|, and reflect the fact that at these points fl^, 
cannot be used to measure time. 

Dark energy density pc{flc)'- The perfect fluid equation 
of motion for dark energy reads; 



Pdo = SH [I + Wc(^c)] Po • 



(11) 



Using Eqs. ^ and ^ in the above equation and drop- 
ping radiation we have that: 



"^dpe 
Po Pc 



dflr 



n,{i~n. 



(12) 

The integrals are the same as those that were required 
for the a(ric) solution, and we find: 



Pci^c) = Po 



rtr 



l-f^p 



Wci^c) 



.Fin.) 



Hubble rate, iJ(J7o): The Hubble rate is simply 

I Pci^c) 



H{n^) = 



3n^ 



(13) 
(14) 



with po(^c) given by equation Eq. (13) 



Angular diameter distance, dA{^c)'- In a flat universe, 
the angular diameter distance to an object at redshift 

z ~ (1/a — 1), is given by: 



dA{z) ^ {I + Z)-' - 

Jt a 



(15) 



Again neglecting radiation, we can use Eq. (l3|, to sub- 
stitute for di to give 

dAi^c) = a(r2o)Go X 

dflo f^?(l - ^cf{wo + wiil^ + W2nl)^ e^ , 

(16) 



where a = -i(Q;-f5), ^ = -l(/3-|-4), 7 = -§(7-1-6) and 

F^-^F. The constant Go = J^n^^il - noy^°{wo + 

wiflo + W2fl^r°, where ao = g(a + 2), /3o = |(/3-2), and 
7o — g7. This integral can be performed analytically if 
W2 = 0: 

Fi (a + 1; -p,-r,d + 2;fl. 



Wl 



n, 



Wo 



- const. 
(17) 

where Fi{...) is the Appell hypergeometric function, and 



G{n,) 



Gq 

a+l 



n^+'w2 . 



(18) 



Up to a constant, this is the final result for the angu- 
lar diameter distance to first order in the expansion of 

Wo(l^c)- 



III. CONSTRAINTS FROM OBSERVATIONAL 
DATA 

In this section we present constraints on the param- 
eters of our dark energy equation of state parametrisa- 
tion Eq. ^ , by performing a global Monte Carlo Markov 
Chain (MCMC) fit to data. We take the set of base pa- 
rameters 



i9== {n„,h^, Ha, Wn}- 



(19) 



and use a modified version of the CosmoMC code j35] to 
sample from the joint posterior distribution of these pa- 
rameters. 



P(i?|a;) 



£{x\'d)r{'d) 



(20) 



where C{x\'d) is the likelihood of the data x given the 
model parameters 1? and 'P{'d) is the prior probability 
density. Other parameters, such as Qq = ^0(2^ = 0) 
may be derived from this base set. We refer to a single 
realization of i?, as a sample. 

The data we use are the compilation of differential-age 
measurements of the Hubble rate H{z) [31j, the latest 
and most precise (local) estimate of the Hubble constant 
Hq [37], and the Union2 SNIa compilation [5S]. These 
data span a redshift range z = — 1.75. 

Our parametrisation is only valid if fio is a monotonic 
function over the redshift range of interest. Monotonicity 
of flc is broken whenever Wo(^c) — X^n ""^"^"(^o) ~ *-*• 
Since we know that We{z = 0) is well constrained to be 
negative by a variety of different observations [IHSl [39] , 
we impose the hard prior wd^c) < at all redshifts of 
interest, z — — 1.75. This corresponds to (l^ > 0. So, 
for a given sample 1?, if any single data point within this 



condition 



monotonic interval 



W2 > 0, fJ- < r^o < f^+ [max(ri_,0),min(ri+,l) 
W2 < 0, n+ < Q.0 [max(r2+, 0), 1] 

iO2<0, f7_>no [0,min(n_,l)] 



U)2 = 0, lUl > 

■U)2 = 0, «;i < 



[0,min(f7_,l)] 
[max(t7_,0),l] 



TABLE I. The intervals where f^c > (if radiation can be 
neglected) for Wciyic) = wo + Wif2e + ^2^6- 



redshift range yields Wc{^c) > then the entire sam- 
ple is rejected and does not feature in the evaluation of 
the likelihood function. For monotonicity to be broken 
within the interval 2; = — 1.75, a fairly rapidly varying 
equation of state would be required, which is disfavoured 
given existing constraints [14]. 

If we truncate the expansion oiwci^c) at second order, 
we can take advantage of our analytic solutions of Sec- 
tioniniwhen numerically implementing this prior. Firstly, 
the roots of Wcl^c) = w^ + wiJ^o + W2^^ are computed: 



0+ = 



— Wl ± \Jw\ — AwqW2 

2wo 



(21) 



We remind the reader that the Wn are related to w„ 
through Eq. (A7). The prior on Wc today becomes 
Wq + wiHq + W2^o < 0. If wf > AwqWi, then depend- 
ing on the values of wq, wi and W2, there are three re- 
gions of monotonic 51c that can be defined. These are 
summarised (along with the case of W2 = 0) in Table p] 
Alternatively, if wf < AwqWi then there are no real roots, 
and We{z = 0) < is sufficient to guarantee that flc re- 
mains monotonic at all times. 

We would like to convert the redshift Zi (or scale fac- 
tor tti) of each data point to a density r^o;. Since Eq. Q 
cannot be inverted analytically, this must be done nu- 
merically. We use a simple bisection method, where the 
boundaries of the search interval over 51c correspond to 
the limits of the monotonic region of 51c given in Ta- 
ble [l] Since 51o is by definition monotonic in this region, 
there will only be one single root, 51o, to the equation 
a(51ej — flj = 0. If the bisection fails to find a solution 
in this interval, then the solution does not exist in this 
interval, and must exist where Wc > 0. If this is the case, 
the sample is rejected, and does not feature in the eval- 
uation of the likelihood function. Those samples which 
generate Wc < for z < 1.75 but Wc > for z > 1.75 are 
still accepted, since our demand is that 51c need only be 
monotonic over the region where our data lie, z < 1.75. 

We emphasise that the mapping from a^ to Slg. does 
not rely on having an analytic solution a(51c). The back- 
ground equations can be easily solved numerically, and 
so the mapping from a; to f2ci is trivial once the re- 
gions of monotonic 51c are known. Hence, in principle 
the Qi — > flci mapping can be performed for an arbitrary 
number of terms in the Chebyshev expansion of Wq. 

We have found the data we use is insufficient to place 



constraints on W2, and so we set it to zero, and constrain 
only the first two terms in the expansion, and so specialise 
to the case: 



i9 = {rtmh^, Ho, Wo, wi} . 



(22) 



Now we summarise the data that we use in our MCMC 
analysis. 

HST Hq prior: The current best (local) measurement 
of the Hubble constant comes from the observation of 
Type la supernovae (SNe la) via the Wide Field Camera 
3 (WFC3) on the Hubble Space Telescope (HST). This 
estimate is Hq = (73.8 ± 2.4)kms"^ Mpc"^ which 
includes systematic errors, corresponding to a 3.3% 
uncertainty [37r\ 

H{z) from differential-age techniques: A weakness of 
supernova observations, BAO angular clustering, weak 
lensing, and cluster-based measurements is that they rely 
on an integral of the expansion history, rather than the 
expansion history itself. The diiferential-age technique 
circumvents this limitation by measuring the integrand 
Az/dtf^ directly, or in other words, the change in the 
age of the universe as a function of redshift. This can 
be achieved by measuring the ages of passively evolving 
galaxies with respect to a fiducial model, and so does 
not rely on computing absolute ages. We use the compi- 
lation of eighteen measurements of Hubble rate Hohs{z) 
that are quoted in |36j . which span the redshift range 
z = — 1.75. For each of the eighteen data points i, 
we convert from redshift Zi, to dark energy density ^^ 
via the bisection algorithm discussed above, and fit our 
samples by minimising: 



Xh{z) 



E 



[g(51e.,1?)-gobsffioj]^ 



(23) 



where Ui are the measurement variances. We use 



Eq. (14) to compute i/(51ci,i?). 



Union2 SNIa sample: We use the Union2 SNIa com- 
pilation released by the Supernova Cosmology Project 
(SCP) [38J, which consists of 557 data points, spanning 
a redshift range z — Q — 1.4. The statistical analysis of 
such SN samples rests on the definition of the distance 
modulus: 

^(z,) = 51ogio[(l + z^fdAiz^, i9)] +25 + ^io, (24) 

where the angular diameter distance (iyi(zi,i9) was de- 
fined in Eq. (15 1. The nuisance parameter /xq, encodes 



^ This local measurement of Hq is in tension with the recent Planck 
measurement [26], Hq = (67.3 ± 1.2) kms"! Mpc"!, from CMB 
data alone. For a discussion of the differences between these 
measurements see |26| . 



the value of Hq^ over which we analytically marginalise 
with a flat prior. This is the standard marginalisation 
procedure, see for example [351 EO] , and is equivalent to 
marginalising over SNIa absolute magnitude. We convert 
from Zi to rioi , and minimise the following expression: 



XSNIa 



T^-1 J d (Cgj^j^) d 



d^Csiiiad 



-1-1 

■"SNIa 



(25) 



The column vector d contains the theoretical minus ob- 
served distance moduli: 



d, = 51ogio[dA(f^o, , "&)a ^(f^c, , i9)] + 25 - A*obs(f^cJ , 

(26) 
whilst CsNia — C'sNia(^ci, ^G ) IS the covariance matrix. 
We use Eq. ([8| to compute a~^(rioii '•^) and perform the 
integral for d^(f2ei , "&) in Eq. \1Q\ numerically (with W2 = 
0). The total y^ to be minimised is then 



2 _ 2 I 2 I 2 

X — XHo + Xh{z) + XsNIa 



(27) 



Taking the base set of parameters (22), we compare 



the performance of our wd^c) parametrisation 

Woi^c) ^ Wo + Wifle , (28) 

against two common parametrisations that may be found 
in the literature: the Chevallier-Polarski-Linder (CPL) 
parametrisation |15l 116) , which is currently favoured by 
the WMAP team [5]: 



.r^^)=-?""+-r" 



(29) 



and the parametrisation of Gerke and Efstathiou 
(GE) 02,: 



,.GE 



(^) 



,.GE 



„GE 



In 



l + z 



(30) 



We begin by assuming flat priors on all base parameters, 
ensuring that they are wide enough such that they do 
not affect the posterior distributions of the parameters. 
We also impose the hard prior Wc{z) < for z < 1.75 for 
the CPL and GE parametrisations in order to facilitate 
a fair comparison to our parametrisation which requires 
this prior in order to be valid. The maximum likelihood 
values of the base parameters and their I — a deviations 
are summarised in the top third of Table III] 

Since we allow Jig (the dark energy density today) to 
vary, our equation of state parametrisation, Eq. (28), ac- 
tually has three free parameters, Wq , Wi and Oq. This is 
different to CPL and GE, which depend on only wq and 
wi. In order to directly compare the three parametrisa- 
tions then, we also show constraints on the derivatives 
oi Wf, at z — 0. For the CPL parametrisation we have 



w, 



CPL I 



w, 



CPL and dw;CPL/dz| 



Wi^^, whilst 



for the parametrisation of Gerke and Efstathiou, we have 
w^^\z=o = w^^ and dw^^/dz\^^o = -wp^. The equiv- 
alent expansion parameters in our parametrisation are 
given by: 



We\z=0 = Wq +^1^0, 

dw, 



dz 



z=0 



3wi{wo + WirJo)^o(l — ^o) 



(31) 



From these two equations, it is clear how our parametri- 
sation depends upon ^Iq. This comparison in terms of 
derivatives of Wc is necessary if we wish to compare like- 
for-like expansion parameters. 

Now, since we have assumed flat priors on all of our 
base parameters, the derived parameters w^^'"''^^|z=o 
and dw^^^''^^ /dz\z^o will also have flat priors, since 
these derived parameters are simply linear combinations 
of the base parameters. For our dark energy cosmic clock 
parametrisation however, the analogous derived parame- 
ters, Eq. (31 ), will not have flat priors since they are non- 
linear combinations of the base parameters. In order to 
make a fair comparison, we adjust the sample likelihoods 
and weights in our MCMC chains for the dark energy 
clock parametrisation in order to obtain flat priors on 



Eq. (31) 



-InC 



-ln£+ln|J| 



weight — > weight X |J| . (32) 



Here |J| = \det{J)\ = | 3{wo + wiflDil -Hq)]. The con- 
straints on Wo 1 2=0 and w^j^^o with both flat and non-flat 
priors are shown in the middle third of Table [III ^^ can be 
seen from the Table, the constraints on 'Wc\z=o are robust 
to changing the prior, whilst w'g\z=o shows a weak depen- 
dence on the prior. To compare the performance of the 
three parametrisations, one must compare constraints in 
the shaded region of Table [ll[ This region displays the 
main results of our paper, which demonstrates that our 
parametrisation performs better than CPL and GE. 

To conclude the discussion of our Table, we summarise 
our results by presenting two Figure of Merit (FoM) cri- 
teria, in addition to the x^ P^r degree of freedom. These 
may be found in the lower third of the Table. The FoM 
criteria are useful statistics that encapsulate the con- 
straining power of cosmological data, given a particu- 
lar dark energy parametrisation. The FoM for the CPL 
parametrisation has been defined by the Dark Energy 
Task Force (DETF) [HI |32] as the inverse of the area of 
the 95% CL region in the wq — wi plane: 



FoMdetf 



1 



v/det C('u;eU=o , <U=o) 



(33) 



where C is the 2D covariance matrix for wdz^o and 
w[.\z=o- Fig. [2] shows these regions for our dark en- 
ergy clock parametrisation with flat priors on We\z=o and 
w'Jz=o (red) compared against the CPL (blue) and GE 
parametrisations (green). Clearly, the better preforming 
parametrisation will have a higher FoMdetFj which we 
find to be our dark energy clock. The second FoM that 





Dark Energy Clock 


CPL 


GE 




flat priors on toq and w\ 


flat priors on TOc|z=o and TOo1z=o 






Ho 

Wo 
lOl 


0.127tro3 
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-1.53tl:f3 

0.74t|S 
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l-78_2'.38 


Oo 
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0.76l°:g? 

— i.iO_Q 24 

0.37lli 


0.671°:°^ 

1 r,i+0.25 
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-0.95tl:?°9 


0.67412:°^, 

-1.14tSi^ 
1 64+2-'76 


0.660lHi 

1 14+0-26 
— -L--L4-0.27 

-1.78t^:fs 


Xdot 

FoMdetf 

FoMwang 


0.945 
2.797 
56.319 


0.936 
2.423 
62.327 


0.945 
1.342 
30.566 


0.945 
1.581 
41.113 



TABLE II. Maximum likelihood parameter values and their 68% CL upper and lower limits. The top third of the table gives 
constraints on the base parameters, whilst the middle third gives constraints on the derived parameters. The lower third gives 
the x^ per degree of freedom, and two FoM criteria, which are useful statistics that encapsulate the constraining power of 
cosmological data, given a particular dark energy parametrisation (see the main text for a detailed explanation of these FoM 
statistics). For our Dark Energy Clock parametrisation, we show constraints on the expansion parameters assuming flat priors 
on Wo and w\ (so non-flat priors on W(,|z=o and 111012=0) and constraints with flat priors on KJe|z=o and Wclz^o- The shaded 
region indicates the parameters which may be directly compared between three parametrisations. 



we quote is: 

FoMwang 



^dct C(wcU=o , WcU=o , ^0) 



(34) 



which was defined in [33]. This extends the two dimen- 
sional FoMdetf to include all the dark energy parame- 
ters. We include VLq = Vtc{z — 0) since in our dark en- 
ergy clock parametrisation, the derivatives d'^Wc/(\-z^\z=Q 
depend rip- As explained in [l^, this FoM penalises 
parametrisations that yield very correlated estimates for 
the dark energy parameters. With flat priors on wq and 
wi (so non-flat priors on We\z=o and w'J\z=o), we find 
that ^6 12=0 in our DE-clock parametrisation is highly 
correlated with i7o as is shown in the left panel of Fig. |3] 
The right panel of Fig. [3] shows the same contours, but 
this time with fiat priors on Wc|z=o and 'w'^\z=o^ the ef- 
fect of which is to reduce the correlation between Wo|z=o 
and ^^12=0- This is reflected in the higher FoMwang in 
Table ini 

As a final way of comparing the three parametrisa- 
tions, we plot in Fig. |4] the redshift evolution of the 
dark energy equation of state for each parametrisation. 
The shaded regions correspond to the values of wdz) 
that are ruled out at 68% CL, and were computed di- 
rectly from the MCMC chains, and not through Gaus- 
sian error propagation. That is, for each sample '& in the 
chain we generate the distributions We{z)^^^ (Eq. 29) 
and u'c(z)^^ (Eq. 30 1 over the redshift range where the 



data lie, z = — 1.75. From these distributions, one 
can easily compute the confidence regions. To generate 
the equivalent distributions using our dark energy clock 
parametrisation, Wc(i^c) = wo+wir^c, we take each sam- 
ple I? and convert from redshift to Vt^ using the bisection 
algorithm described in Section. |TT| 



As can been seen from Fig. |4J the upper 68% CL limit 
on Wc quickly jumps to zero for all three parametrisa- 
tions. This is because the hard prior Wc{z) < cuts off 
the distribution for Wc{z) beyond a given redshift, and 
so the upper 68% CL limit on Wc at these redshifts is 
simply Wc < 0. Whilst the CPL and GE parametrisa- 
tions allow Wo{z) to stray well below —1, our dark energy 
clock constrains Wc{z) to be close to —1 across the en- 
tire redshift range of interest. This is due to two effects. 
The first is the lower value of w'g\z=Oj (or more precisely 
wq) for our parametrisation compared to CPL and GE. 
The second, more important effect is that our expansion 
parameter ^dz) decays much faster with increasing red- 
shift compared to the expansion parameters z/{l + z) 
and In [1/(1 -I- z)] of the CPL and GE parametrisations 
respectively. For example, we find ^o{z = 1) « 0.05, 
whilst z/{l + z) = 0.5 and In [1/(1 + z)] = -0.69 at 
z = 1. Both of these effects keep Wc{z) close to —1. 
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FIG. 2. 2D 68% (dark shading) and 95% (light shading) marginahsed contours in the Wc\z=o-w'c\z=o plane. Left panel: Red 
regions: DE-clock with flat priors on uio|z=o and w'J\z^o\ blue regions: CPL parametrisation. Right panel: Red regions: 
DE-clock with flat priors on Wc\z=o and u;^|z=o; green regions: GE parametrisation. 
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FIG. 3. 2D 68% (dark shading) and 95% (light shading) marginalised contours in the u;e|z=o-f^o plane. In both panels red 
regions correspond to the DE-clock and blue regions to CPL. Left panel: Flat priors on wo and lui. Notice that 'Wc\z^o is 
highly correlated with r^o. Right panel: Flat priors on Wojz^o and w'J\z=o. The correlation is reduced. 



IV. COMPARING TO THEORY 

As discussed in Section |ll| since we expand in a ba- 
sis of orthogonal Chebyshev polynomials, any dark en- 
ergy model which predicts a monotonic equation of state 
'Wtheory(^c), Can be dircctly related to our parametrisa- 
tion. We illustrate this by placing constraints on the 
Chebyshev expansion coefficients (the w„ of Eq. (l5])) of 
two popular quintessence models, the single exponential 
potential (lEXP) [291 [30] and the supergravity inspired 
potential (SUGRA) [SJ. Using the 7?o, H{z) and SNIa 
data discussed in Section |III[ we vary the base set of pa- 
rameters 



1? 



theory 



= {Q^mh , H{), A} , 



(35) 



of both quintessence models, and use a modified version 
of the CosmoMC code [35] to sample from the joint pos- 



terior distribution of these parameters. For each sam- 
ple i^thoory in the MCMC chains, we can construct the 
equations of state, wiexp(^c) and wsuGRA(f^c)- Once 
these functions are known, the expansion coefficients, 
Wn, of each sample may be extracted by appealing to 
Eq. ([5|. If a sufficient number of samples are taken, 
we can generate distributions for the Wn- The condi- 
tion that Vte{t) must be a monotonic function of time, 
results in the prior M;iEXp(^e) and ■u;suGRA(^e) < over 
the redshift range of interest. We choose the region over 
which the Chebyshev polynomials are orthogonal to be 
Q,c = [^c{z = 1-75), r^o], which spans the redshift range 
of the data. 



In Table III we quote the maximum likelihood values 
and the marginalised 68% confidence limits of the 
Chebyshev expansion coefficients w„ for the lEXP 
and SUGRA models. We also give constraints on the 






FIG. 4. The redshift evolution of the dark energy equation of state Wc for dark energy clock (red), CPL (blue) and GE (green) 
parametrisations. The shaded regions correspond to the values of Wo(z) that are ruled out at 68% CL. Notice that the upper 
68% CL limit on w^ quickly jumps to zero for all three parametrisations. This is because the hard prior Wc{z) < cuts off the 
distribution for Wc{z) beyond a given redshift, and so the upper 68% CL limit on We at these redshifts is simply Wc < 0. Whilst 
the CPL and GE parametrisations allow We{z) to stray well below —1, our dark energy clock constrains We{z) to be close to 
— 1 across the entire redshift range of interest. 





lEXP 


SUGRA 


Wo 


-0.975 (-1.0,-0.972) 


-0.726 (-0.756,-0.696) 


Wi 


0.01 (0.0,0.012) 


-0.093 (-0.104,-0.092) 


Wo 


-1.006 (-1.008,-1.0) 


-0.298 (-0.312,-0.20) 


Wi 


0.069 (0.0,0.080) 


-0.770 (-0.90,-0.758) 


Wc\z^O 


-0.954 (-1.0,-0.947) 


-0.912 (-0.930,-0.893) 


"'cU^O 


-0.03 (-0.04,0.0) 


0.339 (0.298,0.379) 



TABLE IIL Maximum likelihood values and marginalised 
68% confidence limits (lower , upper) for the first two Cheby- 
shev expansion coefficients corresponding to the lEXP and 
SUGRA quintessence models. We also quote constraints on 
the derived parameters w„ (which are related to the w„ by 
Eq. (A7l), and the equation of state today, Wc\z=o, and its 
derivative, w'Jz=o- 



derived parameters Wc\z=o and w[.\z=o- We compare 
these constraints in Fig. [5] by superimposing the 2D 
marginalised contours in the 'Wc\z=o~Wc\z=o plane for 
the two quintessence models upon the contours of our 
dark energy cosmic clock. 



As can be seen from the figure, both quintessence models 
are consistent with the constraints on the dark energy 
equation of state at 68% CL using our parametrisation. It 
is in this way that models of dark energy can be compared 
to our parametrisation, much like theoretical predictions 
of selected inflationary models superimposed upon the 
marginalized confidence regions for n^ and ro.oo2 using 
the recent Planck data ^E\. 

We note that even though we have assumed flat pri- 
ors on all base parameters, the priors on the expansion 



coefficients Wn will not be fiat. This is because the in- 
tegral in Eq. (I5| will in general be some complicated 
function of the base parameters. Furthermore, the lim- 
its QT^ = n^{z = 1.75) and n^""" = ^o of the integral 
Eq. ([5| will themselves have probability distributions. A 
more careful and detailed comparison involving identi- 
cal priors will be left to future work, but we include this 
simpler study here in order to give the reader a general 
impression. 



V. DISCUSSION AND CONCLUSIONS 

Various astrophysical constraints suggest that the dark 
energy density f2c niay have increased monotonically 
throughout cosmic history [21 |25l |44j. Acknowledging 
this apparent monotonicity we have introduced a com- 
pletely new parametrisation of the dark energy equation 
of state which uses the dark energy density ftc{t) as a cos- 
mic clock. Our parametrisation has several advantages, 
perhaps the most important being that flc is a physi- 
cal quantity, directly related to the properties of dark 
energy. Furthermore, assuming < fie < 1, it makes 
for an ideal expansion parameter since f2c is a naturally 
small number. By fitting to SNIa and H(z) data, we have 
demonstrated that our parametrisation performs better 
than the popular Chevallier-Polarski-Linder (CPL) and 
Gerke and Efstathiou (GE) parametrisations, and that it 
is robust to the choice of prior. Expanding in orthogo- 
nal polynomials also allows us to relate models of dark 
energy directly to our parametrisation, which we have 
illustrated by placing constraints on the expansion coef- 
ficients extracted from two popular quintessence models. 

The dark energy density f^e can only be used as a cos- 
mic clock if it is a monotonic function of time. As can be 
seen from Eq. ([3]), for a universe containing dark energy, 
pressureless matter (tUm = 0), and radiation {w.y = o), 
there are turning points in il^ whenever w^ = ^ m the 
radiation dominated era, and in the matter era when- 
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FIG. 5. The 2D 68% (dark shading) and 95% (hght shad- 
ing) marginalised contours in the Wc\z=o-w'c\z^o plane for 
the lEXP and SUGRA quintessence models superimposed 
upon the corresponding contours of the dark energy clock 
parametrisation. 



ever Wc — 0. Many scalar field models of dark energy 
possess scaling solutions on which the dark energy equa- 
tion of state Wc 'tracks' that of the dominant background 
component (see e.g. [29l [30]). In such models, mono- 
tonicity of fie would be spoiled as Wg evolves through 
the radiation and matter dominated eras and transitions 
towards Wo '^ — 1 today. In this paper, we have used 
our parametrisation to probe the dynamics of dark en- 
ergy at low redshifts, z = — 1.75. It is presumably safe 
to assume that J7e remains monotonic over this redshift 
range, since to break monotonicity between z = and 
z ~ 2 would require a rapidly varying Wg-, which is not 
favoured by existing analysis [T31[53]. Hence, in our fit- 
ting to data we imposed the hard prior Wo{z) < for 
z = 0-1.75. 

To probe the high redshift behaviour of dark energy 
this prior would need to be removed, since it would be un- 
realistic to say with any sort of certainty that Wc{z) < 
throughout the entire cosmic history. Hence, to accom- 
modate high redshift data the non-monotonicity of J7c 
would need to be accounted for. This could be achieved 
by piece-wise parametrising Wc in regions of monotonic 
r^c- For example, if radiation can be neglected, these re- 
gions are defined by the roots of the polynomial wd^c) — 
J2'^ w„C/„(51c) = 0. If Wo(rie) = wo + wiflc and wi > 0, 
then there would be two distinct regions: r2o > for 
r^o < —wq/wi (I) and f2o < for il^ > —Wo/wi (II). One 
would then write 



Wci^c) = 



W, 



(11) 









< 






(36) 



and so there would be four free parameters in total. To 
accommodate CMB data, radiation can not be neglected, 
but such regions of monotonicity can still be defined: one 
would need to compute the roots of Eq. ^ exactly, which 



could be performed numerically. 

Finally, it is interesting to make the connection be- 
tween our parametrisation and the dark energy 'flow pa- 
rameter' 

F EE K^^ , (37) 



We 



that was introduced in [351 [3^ (see also [3Z])- Here, 
A = —V^/V, where 1^(0) is the scalar field dark energy 
potential and V^ is the derivative of V with respect to 
(/). By considering general dark energy models where the 
field either accelerates or decelerates down its potential 
toward its minimum (dubbed 'thawing' or 'freezing' field 
evolution [48]), the authors of [45] were able to demon- 
strate that F remains nearly conserved until quite recent 
times, z ~ 1 — 2, after which dark energy finally begins 
to take over. This is despite Wo, flc and A all being dy- 
namical. This constant nature of the flow parameter is a 
direct consequence of the fact that the dark energy fleld 
does not exist in a vacuum: instead it has been influenced 
by the long periods of radiation and matter dominated 
epochs prior to the current day. 

If the parameter A is a constant, (as is the case for ex- 
ponential potentials) or remains approximately constant, 
then we have Wc ^ —1 + FX^flc, which looks very much 
like our dark energy clock parametrisation with wq = — 1 
and Wi = FX^. This indicates that, so long as F re- 
mains approximately constant throughout the long radi- 
ation and matter dominated eras, then the behaviour of 
a wide range of scalar field dark models should be well 
captured by We(^c) = wq + wiflc. 
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Appendix A: The Chebyshev Polynomials of the 
Second Kind 

The Chebyshev polynomials of the second kind are de- 
fined by the recurrence relation 

C/o(x) = l, Ui{x)^2x, 

U„+i{x) = 2xUnix) ~ U^^iix) , (Al) 

and obey the following orthogonality condition 

Un{x)U^{x)Vl - x'^dx = -(5„™ . (A2) 



11 



We can shift the interval over which the polynomials are 
orthogonal by choosing a suitable inner product. Let 
Unix) = Un{f{x)) where f{x) is monotonic on [a, 5], and 
satisfies /(a) = —1 and f{h) = 1. The simplest choice is 



/(^) = 



2x ~ a ^ b 



(A3) 



Then, the orthogonahty condition becomes 

Un{x)Um{x)\/{x — a){b — x)dx — —{b — a)^S„ 



(A4) 

Any function g(x) which is continuous in the interval 
of orthogonality [a,b], may be expanded as a series of 
Chebyshev polynomials: 



g{x) = y^^w„U„{x) , 



(A5) 



From the orthogonality condition, Eq. (A4) we have 



Ab-a^Ja 



g{x)Un{x)\/{x — a)[b — x) dx . 



(A6) 

Under certain conditions of the interpolated function 
g{x) (Dini-Lipschitz continuity), the Chebyshev interpo- 
lation converges when the number of nodes tends to in- 
finity. For numerical implementation, it is convenient to 



make the change of variable x 
16 r-n 



(b — a) sin^6' giving 



Wri 



g [a + {b ~ a) sir 



X Un (a + (fe — a) sir 



sin^^cos'^^ 



de. 



In the case where q{x) = wd^c), and b = 51™^'^, a = 
r2""", we can write Eq. ( A5 1 at second order as Wo{^c) — 
wq + wifle + W2^e- The Wn are given in terms of the Wn 



as: 



Wo 



Wo 



2{n- 



f^max „ f^min 
in ^ f7max)2 



f^max-) 

Wi 



1 



Wi 



W2 



(^max _ f]min)2 

4 . 16(f7, 

^max _ f^min ^1 ( J^m 

16 

;W2 ■ 



W2 






i^^e 



^miny 



(A7) 



Appendix B: The Sic— clock particular solutions 

The powers a, /3 and 7 that appear in the solutions for 
a(51c) and Pe(^e) (Eqs. ((sl) and (l3^ ) are combinations 
of Wq, Wi and W2, and correspond to coefficients of the 
partial fraction expansions of Eqs. (It]) and (12). As a 
result, if any one or more of the coefficients wq, wi or 
W2 is exactly zero, then the powers a, /3 and 7 and the 
function F(J7o) will change. In table IV we list all the 
possible solutions for different combinations of wq, wi 
and W2, where one or more of them is zero. 
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TABLE IV. The values of the powers a, /3 and 7 and the function F(f2c) for particular cases of zero wo, w\ or u)2. Here, 
g — \J \wq'W2 — wf, and wt = wo + wi + W2. Notice that F{Qe) may always be kept real by using the relation arctan(i2:) = 
iarctanh(2:) when necessary (for example if 4iuo'U'2 < wl in the definition of q). 
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